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Abstract 

Motivated by a problem arising from pharmaceutical science [B. Baeumer et al, Discr. Con- 
tin. Dyn. Sys. B 12], we study random walks on the contact graph of a bidisperse random sphere 
packing. For a random walk on the unweighted graph that terminates in a specified target set, 
we compare the number of steps and the euclidean distance covered. We find a linear relationship 
between the two metrics with a proportionality constant that can be calculated from the edge 
length probabilities of the contact graph. 
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I. INTRODUCTION 



Recently, we developed a model to predict the release kinetics of matrix tablets [1]. A 
matrix tablet is a device used to deliver a drug and to release it in a controlled fashion over an 
extended period of time. They are formulated as powder mixtures of the pharmaceutically 
active drug and at least two inactive ingredients, namely water soluble excipients and water 
insoluble polymers. Upon placement of the compressed tablet in fluid, the polymer matrix 
remains largely intact while soluble excipient and drug particles are dissolved and carried 
away by diffusion. It is of great interest to predict the time course of the drug release and its 
dependence on the composition of the powder mixture. In [1] we proposed a mathematical 
model for the drug release process, based on a discrete random walk model which serves as 
the starting point for the present paper. 

In [1], we begin with the creation of a random dense sphere packing V, see [2, 3, 7, 8] 
for recent discussions of this notion. Here, it will be the result of the Lubachevsky-Stillinger 
(LS) protocol [4, 5]. Briefly, growing spheres move and collide until the packing fraction 
increases only marginally between two collisions and the procedure is halted (alternatively, 
the pressure, the sum of the squared velocities exceeds a given threshold). From the random 
sphere packing so obtained, we construct the contact graph by putting an edge between 
sphere centers that are within a certain distance of each other. In our implementation in 
[1] we label some particles as "drug" particles and seek their escape routes to "exterior" 
particles. We count the number of steps that it takes each particle to reach one of the 
exterior vertices, repeating the random walk for each starting point a certain number of 
times. From the histograms of the number of steps we predict the cumulative release profile 
of the simulated matrix tablet, in good qualitative agreement with experimentally formulated 
matr.ix tablets As we considered only monodisperse packings in [1], no distinction between 
topological and geometrical random walks was necessary. However, it is well possible that 
the powder particles of different compounds have different sizes. This leads to the question 
whether there is a relationship between the two metrics. 
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II. TOPOLOGICAL AND GEOMETRICAL RANDOM WALKS 



Let V C [0, l] 3 be the set of centers of a bidisperse random sphere packing on the flat 
torus, where the spheres have radii r\ < r 2 with m 1; respectively m 2 spheres in each class. 
Let vertices Xj and Xj be the centers of spheres with radii r l and r J , respectively and fix e > 0. 
Then Xj and x^- are joined by an edge in the graph G £ = (V, E e ) if |xj — Xj\ < (1 + e)(r t + r J ). 
The infimum of all e > such that G £ is connected, is called the level of connectedness. 

Let V a and be subsets of V, chosen randomly among all vertices V. We perform random 
walks on the graph G £ = (V, E e ) that start in V a and terminate in V w . For simplicity, we will 
work with equal weights for all edges in this paper (but see the Discussion in section IV). 
For such a random walk we count the steps (the "topological" metric) and it's euclidean 
length (i.e. the sum of the lengths of all the edges in the walk, the "geometric" metric). 

Let tti, 7Ti 2 and 7r 2 be the probabilities for an edge e e E e of G e — (V, E £ ) to be (approx- 
imately) of length 2r 1; r± + r 2 and 2r 2 , respectively. Then we expect a linear relationship 

l(n) x (2ri7Ti + (n + r 2 )7Ti2 + 2r 2 7r 2 )n =: a th n, (1) 

in the sense that the quotient of the two expressions approaches 1 as n — ► oo. 

We use an implementation of the Lubachevsky-Stillinger protocol developed by the Com- 
plex Materials Theory Group at Princeton University [6] with periodic boundary conditions. 
Unless stated otherwise, we fix the total number of spheres to be 400. For the volume frac- 
tion p of small sphere volume vs. the total sphere volume we have p = 1+ g CT3 , where £ = ^ 

and a = —. 

n 

III. RESULTS 

We consider the correlation between the packing fraction <ft and the level of connectedness 
of a packing. It is very hard to pack tightly mixtures of a small aspect ratio a -1 = ri/r 2 
and a small volume fraction p of small spheres. Geometrically, an individual small sphere 
may float freely in the interstitial space between the large spheres, even though the overall 
packing fraction may be significantly higher than the packing fraction of a monodisperse 
packing (see figure 1, left panel, "+"). 

After creating the sphere packings, we select those that have a level of connectedness 
< 0.05 for the further investigations. In figure 1, right panel, we show the probabilities 
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FIG. 1: (Left) Plot of the level of connectedness vs. the packing fraction of 50 respectively 
100 packings for two choices of aspect ratios cr _1 and volume fractions p. The termination criteria 
used in the Lubachevsky-Stillinger protocol are identical. An extreme case of high packing fraction 
but poor level of connectedness is indicated by an arrow. (Right) Probabilities for an edge to be 
short (wi, solid lines) respectively long (-7T2, dashed lines) as functions of the volume fraction p of 
small spheres. 

7Ti and 7T2 of short and long edges, respectively, as functions of o~ _1 and p. 

For each packing we create 10 4 random walks starting at a random initial point and 
terminating in a random set of varying sizes. Each walk passes through a certain number of 
short {2ri), medium (r x + r 2 ) and long (2r 2 ) edges. These numbers are recorded to compare 
them with the theoretical edge length probabilities tti, 7ti2 and 7T2 that are used in equation 
(1). We select all walks that share the same number of steps n and take the average of their 
euclidean lengths, this gives the quantity l(n). When we plot this against n, we observe a 
linear relationship with slope a exp for small values of n, before the points become scattered, 
see figure 2. The theoretical value a t h from (1) overestimates the actual slope somewhat, 
but the ratio a^/a exp between theoretical and experimental slope decreases with the size of 
the target set (and hence an increased average duration of the random walks), see figure 3. 
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FIG. 2: (Left) Plot of the average lengths of random walks of n steps on a packing of 400 spheres 
with cr _1 = 1/4 and a volume fraction p = 33%. A total of 10 4 random walks were created, ending 
in a target set of size 50. The solid line shows a linear fit to the data (giving the slope a exp ), 
while the dashed line has the slope predicted by equation (1), where the probabilities tti, tt\2 and 
7T2 are taken from the contact graph. The dashed-dotted line uses the same equation, but with 
the empirical probabilities taken from the random walks. (Right) As in the left panel, but with 
\V U \ = 10. 

IV. DISCUSSION AND CONCLUSION 

The random walks considered here mimic the diffusion of a drug molecule from its original 
position to the edge of a matrix tablet [1]. The length distribution of these walks is used 
to calculate the escape times for molecules and the cumulative distribution function of the 
escape times is taken a prediction for the release profile of the tablet. If the particles of the 
powder mixture have different sizes, the definition of the "length" needs to be revisited. 

Here we have provided evidence that the geometric and topological metrics on the contact 
graph of a random dense sphere packing are related by the formula (1), in the limit of long 
walks that sample correctly the edge length distribution. While the average ratio a t h/a exp 
for the walks with size of the target set |V^| = 10 is 1.16 (figure 3, "+"), we expect this 
ratio to approach 1 as the number of spheres in the packing increases. This is the focus of 
ongoing research. 

The probabilities n\ and 7T2 are determined from contact graphs for different types of 
sphere mixtures. As figure 1 (right panel) suggests, these functions are monotone in both 
the aspect ratio a~ x and the volume fraction p of the sphere mixture. We conjecture the 
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FIG. 3: A plot of the ratio ath/ a exp of theoretical and experimental slopes for random walks on 93 
packings, with different sizes of the target sets. 

following limits 

if p = J 1 if p = 
lim vr 2 (a ,p)= < 

1 if p > 0, ff ^^ 0+ [ if p > 0, 

It is difficult to make a prediction for a -1 — > 1~, as then the distinction between short, 
long and medium edges disappears. Notice the change from concave to convex in 7Ti as a -1 
increases from 1/4 to 3/4 in the right panel of figure 1, while no change of convexity is seen 
in the corresponding values of 7T2- 

A natural generalization is to consider fc-disperse packings with k > r radii classes. 
Then there are k ^ k ^ possible edge lengths in the contact graph and a straightforward 
generalization of equation (1) is expected to hold. The situation is much more complicated 
if the radii of the spheres themselves are a random variable. 

In the present work we have considered only the case of random walks on a graph where 
all edges carry the same weight. The motivating problem from pharmaceutical science in 
[1] however required that some edges are much harder to access than others, namely those 
ending at particles of polymer type. Although the edge length distribution (tti, 7Ti 2 , 7r 2 ) is the 
same, a long random walk is then no longer going to sample all edges with these probabilities. 
The conjectured length of random walks (1) will need to take that into account. 



lim 7ri(<7~ 
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